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Abstract 



V) , We consider the possibility that tunneling in a degenerate double- well poten- 

tial in de Sitter spacetime leads to coherent oscillations of quantum probability 
to find the system in a given well. We concentrate on the case when the mass 
scale of the potential is much larger than the Hubble parameter and present 
a procedure for analytical continuation of the "time" -dependent instantons, 
which allows us to study the subsequent real-time evolution. We find that the 
presence of the de Sitter horizon makes tunneling completely incoherent and 
calculate the decoherence time. We discuss the difference between this case 
and the case of a ^-vacuum, when tunneling in de Sitter spacetime preserves 
, quantum coherence. 



Ph! frep-ph/0111194 

*f; : 1 Introduction 



Tunneling in de Sitter spacetime had been a subject of many papers 0-0. Our 
purpose in revisiting it here is to explore its possible connection with decoherence 
of quantum fields during cosmological evolution. Decoherence refers to the onset of 
classical behavior in a quantum system or, more formally, to vanishing of the off- 
diagonal elements of the reduced density matrix. It is well-known that the flow of 
quantum modes of a scalar field through a de Sitter horizon can result in an effectively 
classical behavior of the field's large-scale, coarse-grained component, a phenomenon 
concisely described in Starobinsky's stochastic approach |4j]. 

In laboratory systems considered as candidates for macroscopic quantum coher- 
ence (MQC) experiments, decoherence is attributed to interactions with the environ- 
ment. In this case, it is known that a sufficiently strong such interaction can entirely 
wash out any coherent quantum oscillations associated with tunneling [[J . It is by no 
means obvious that the special dynamics of super-Hubble modes in de Sitter space 



can be interpreted as decoherence in the usual MQC sense |10| . Nevertheless, one may 
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reasonably ask what effect the de Sitter horizon has on various specific manifestations 
of quantum coherence. 

One manifestation of quantum coherence, important in particle physics, is the 9- 
vacuum structure |L1| in gauge theories, such as QCD. ^-dependence of the partition 
sum is a result of interference between topologically distinct paths in the configuration 
space and has no counterpart in classical theory. In theories that are conformally 
invariant at the classical level, such as massless gauge theories in four dimensions, it 
is easy to construct the corresponding instanton solutions for the de Sitter case and 
find a nonvanishing, although — as explained in the conclusion — probably reduced, 9 
dependence. (This implies in particular that the de Sitter expansion alone cannot 
wash out an unwanted 9 angle.) 

Here we want to consider another class of problems where destruction (or per- 
sistence) of quantum coherence is important. These arise in cosmological scenarios 
(such as eternal inflation [|T^]) in which there is an infinite supply of universes, each 
with its own values for at least some of the fundamental constants. The question is if 
there is a mechanism that causes our universe to be in a state with definite values for 
these constants, rather than in a quantum superposition of such states. If one accepts 
that values of the fundamental constants are set by expectation values of some scalar 
fields, e.g. superstring moduli, the question can be reformulated in terms of these 
expectation values. 

In a scalar theory with multiple potential minima, the tunneling rate between the 
minima will be nonzero in a de Sitter spacetime even if it is zero in the flat space. 
An example is provided by a scalar potential with two degenerate minima (although 
many of our conclusions will apply to the nondegenerate case as well). Because the 
Euclidean counterpart of a de Sitter spacetime is compact (a sphere) [[□[], such a 
theory always has finite-action Euclidean solutions — instantons. One such solution 
is the Hawking- Moss (HM) instanton JIJ, for which the field is at a maximum of the 
potential. 

The HM instanton can be used to calculate the rate of transitions between different 
potential minima in the case when the mass of the scalar is smaller than some critical 
value of order of the Hubble parameter H fll], |3[]. This is precisely the case when 
the amplification of super-Hubble modes leads to stochastic dynamics of the coarse- 
grained field ||. The connection between the two approaches has been discussed in 
the literature [|], We have little doubt that in this case the stochastic dynamics 
will cause a rapid decoherence between the minima. 

In the present paper, we concentrate on the opposite case, when the field has 
a degenerate potential but the mass scale of the potential is much larger than H. 
Note that for a degenerate (or nearly degenerate) potential the large mass does not 
imply that we are headed towards the limit |14| of small, nearly flat-space instantons. 
Indeed, we have already noted that for an exactly degenerate potential there are no 
instantons in the flat space. 

In perturbation theory, there is no amplification of super-Hubble modes of a heavy 
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field, and the usual form of stochastic description does not apply. One may wonder if 
in this case tunneling leads to the typical MQC phenomenon of quantum oscillations 
of the probability to observe the system near a given minimum. Our main result is 
that it does not: the presence of the de Sitter horizon makes tunneling completely 
incoherent even for fields that are much heavier than H. 

The plan of the paper is as follows. In Sect. 2, mainly to introduce various 
coordinate systems useful for studying instantons in de Sitter space, we consider a 
case when tunneling does exist in the flat spacetime — the A0 4 theory with the "wrong" 
sign of A. This theory is classically conformally invariant, and the instanton is easy 
to find: it is a conformally transformed Fubini-Lipatov [[L5| instanton. This solution 
has the conventional large- and small-size limits. 

In Sect. 3 we introduce the 4 theory with a degenerate double- well potential 
and discuss the crossover in the properties of the lowest-action instanton solution 
at fi = /i cr ~ H. As explained above, the small-instanton limit does not exist in 
this case, so at fi ^> H, the lowest-action solution is a different brand of "time"- 
dependent instanton: in the "static" coordinates, its "time" duration is order yT x , 
while the spatial extent is of order H^ 1 . In Sects. 2 and 3 we obtain several explicit 
expressions for the turning points of the "time" -dependent solutions. The existence 
of turning points is indeed their key feature, which makes these solutions analogous 
to periodic instantons [|16| in massive gauge theories in flat space (while the HM 
instanton, from this point of view, is analogous to the sphaleron [p"7[ ). 

Sects. 4 and 5 contain the main results of this paper. In Sect. 4 we describe an 
analytical continuation of the "time" -dependent solutions that allows us to consider 
the subsequent real-time evolution. We find that by the time the field that started 
at a turning point reaches a section of constant conformal time rj = i] , it turns into 
a super-Hubble droplet of the new phase (e.g. <fi = ji nucleated inside the phase 
with cf) = —fi). We argue from causality and confirm by numerical integration that 
in de Sitter spacetime such a droplet will expand forever, even if the corresponding 
flat-space theory does not have any expanding droplets (as is the case for a strictly 
degenerate potential). In Sect. 5 we show that, as a result, the tunneling between (f> = 
fi and cf) = —fi is completely incoherent and calculate the corresponding decoherence 
time. Sect. 6 is a conclusion. 

Throughout the paper, we neglect backreaction of instantons on the metric. This 
is a good approximation when the variation of the potential during an instanton 
transition is much smaller than H 2 Mp { . For example, for the double potential, the 
variation is at most of order /x 4 /A, where A <^ 1 is the coupling constant. To consider 
fi ^> H, we then have to require H 2 / A <C Mp 1; a relatively mild restriction. 
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2 Conformally invariant model 



2.1 Instanton solution 

Our first example is the massless A0 4 theory with a "wrong" (negative) sign of A, 
conformally coupled to gravity. The Euclidean action of the field is 



<t> 4 



X 



where is the metric, and R is the scalar curvature of a four-sphere (the Euclidean 
counterpart of the de Sitter space [13|]). The radius of the sphere equals the Hubble 
radius, which in this section we will set to one by a choice of the length unit. Then, 
R = 12. 

There are several convenient choices of metric on a unit four-sphere. One is the 
conformal form 

g^u(x) = - ^ 2 , (2) 

where x is a flat-space four-vector, and x 2 is its length squared. Thus, the sphere is 
related to the flat space by a conformal transformation 

9^(x) -> g^(x) = 7 2 (x)g^(x) (3) 

with 

7 (x) = (l + x 2 )/2. (4) 
The action is invariant under transformations (|3p provided we also transform the field: 

4>{x) — > 4>{x) = •y~ 1 (x)(j)(x) , (5) 

and use the scalar curvature corresponding to the new metric: R — > R. 

The theory in the flat space has a well-known solution, the Fubini-Lipatov instan- 
ton IfE!)] , which describes decay of the false vacuum = 0: 

where b and Xo are arbitrary (real) parameters. Inverting the transformation fl5|), 
with 7 given by ([§), we obtain the corresponding solution on the sphere 



V2b(l + 



x 



2\ 



^ ^ (x — xo) 2 + b 2 ^ 

This solution was easy to find in the conformal coordinates (0), but to study some 
of its properties it is convenient to consider other coordinate systems, which we will 
also use in the following sections. One of these new choices is the system of four 
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spherical angles ((,ip,d,4>), which is most easily introduced if we first go from the 
coordinates x to the corresponding spherical coordinates, by writing the line element 

as 

ds 2 = g^dxV = + 4 ^ 2)2 (dTZ 2 + n 2 dn 2 3 ) , (8) 

where TZ is the length of x, and dO| is the metric on a unit three-sphere: 

dVlj = dip 2 + sin 2 i) d9 2 + sin 2 ip sin 2 9 dcf) 2 . (9) 

We then define the fourth angle ( via 

K = tan(C/2) . (10) 

The line element becomes ds 2 = d( 2 + sin 2 (dQ^. 

The form of solution (|?p in the new coordinates is somewhat complicated, except 
for the case Xq = 0, on which we now concentrate. (This implies no loss of generality 
for our present purposes, as we are free to center the stereographic projection (fj) on 
the instanton.) Then, 

= sin 2 (C/2) cos 2 (C/2) " (U) 

We can see how the solution changes as we change b from b <C 1 to b ^> 1. At small 
b, the solution is concentrated in a small region near ( = (the north pole of the 
four-sphere), where 

This is equivalent to the flat-space instanton (§) with a rescaled size: b — > 2b. As 
b increases, the solution occupies a larger region, until at b = 1 it becomes uniform 
over the sphere. This is, of course, the Hawking-Moss instanton of the model: 

fo = V2- (13) 

Any solution with b > 1 can be obtained from a certain solution with b < 1 by using 
the transformation b — > 1/b, ( — > 7r — (. Thus, at large b the solution is concentrated 
in a small region near the south pole and is approximated by the flat-space instanton 
with b — > 2/6. Therefore, both large and small b limits correspond to instantons 
of sizes much smaller than the curvature radius. Such instantons are practically 
unaffected by the curvature (the Coleman-De Luccia limit [Q). The point 6 = 1 
corresponds to the largest size and strong curvature effects. 
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Figure 1: Coordinates r and r defined by (|T4]), (|T5|). The boundary of the square 
corresponds to r = 0, and the center to r = 1. 

2.2 The static coordinates 

Of particular interest is the form of the solution in the so-called "static" coordinates, 
which are obtained by leaving 9 and <fi of the previous set intact but replacing ( and 
ip with 

r = sin if) sin ( , (14) 
r = arctan(cos-0tan£) . (15) 

The line element now takes the form 

ds 2 = (1 - r 2 )dr 2 + + r 2 dn\ , (16) 

1 — r 2 

where dfl 2 , = d6 2 + sin 2 9 d(f> 2 is the metric of a unit two-sphere. Each of the angles 
( and ip runs from to 7r. For the mapping flU]), ( |T5] ) to be continuous, r should be 
periodic with period 2ir; r runs from to 1. To help visualize this mapping we show 
it pictorially in Fig. [l]. 

Metric (|T6| ) is "static", i.e. independent of r, so it can be easily continued to real 
time. The simplest way to do that is to make r purely imaginary, while leaving r 
untouched. Although this choice of coordinates is not convenient for studying the 
real-time evolution, the possibility to make it shows that turning points with respect 
to r, i.e. configurations satisfying 

<9 r 0(r, r = r ) = (17) 
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for some ro and all r, are of special importance, since they can be used as initial 
conditions for the real-time evolution (cf. ref. [§). 

Solution (|7|) in the static coordinates takes the form 

(p[x) 



{b 2 + 1) + {b 2 - l)VT=r* cosr v ' 

We see that for any b different from 1, the instanton oscillates in the Euclidean time 
t, and there are two turning points, at r = and r = ir. 

When b is close to 1, we can expand fll8|) in small b 2 — 1 = 2e: 



<f>(x) = 0o (l - eVl - r 2 cos r) + 0(e 2 ) . (19) 

So, in this limit, the oscillations are small harmonic oscillations about <ft = <f) . 

Because the model ([I]) is classically conformally invariant, the Euclidean action 
will have the same value for any value of b. Which values dominate the false vacuum 
decay will be decided by quantum correction to the amplitude. (Because the A0 4 
theory with a negative A is asymptotically free, we expect that the quantum correc- 
tions will favor the largest instanton sizes, i.e. b ~ 1.) Certain values of b will be 
favored classically if we deform the model away from the exact conformal invariance, 
by adding a mass term m 2 s (j) 2 /2 inside the brackets in ([!]). Here 

m 2 cS = m 2 + (^-l/Q)R J (20) 

and m 2 and £ are parameters. If \ml s \ <C H 2 , the change in the action can be 
computed using the instanton solutions of the conformally invariant theory, cf. ref. 
fl8| . One can then show that m 2 E > favors instantons with 6^0 and b — ► oo, 
while m 2 s < those with b ~ 1. This is a variant of the well-known crossover [l], [|[ 
in the properties of the dominant classical solutions, and we do not go into details of 
the calculation here. 



3 Double-well potential 

We now turn to models in which, unlike the model of the previous section, there are 
two degenerate or nearly degenerate minima. In the next two sections, we will try to 
answer the main question of this paper, namely, whether in such cases tunneling in 
de Sitter space can lead to coherent oscillations of the probability to find the system 
near a given potential minimum. 

The simplest model of this kind is a scalar 4 theory with an exactly degenerate 
double-well potential. The Euclidean action of the field is 

Se = \J ^^v^{^W0 + ^0 2 -^o0 2 + ^ 4 } , (21) 
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where A is now positive, and the parameter £ can take any value. As before, we 
measure lengths in units of the Hubble radius, so R = 12. The 2 terms in fl2lD can 
be assembled into a single 2 term, with the mass parameter 

/i 2 = fi 2 - 12f . (22) 

We will assume that fi 2 is positive. 

In this section, we will use metric in the form (||), with an additional definition 

y = kiK. (23) 

Thus, the line element is 

ds 2 = — -J- (dy 2 + dQ 2 ) . (24) 
cosh y v ' 

(which is conformally equivalent to the metric of a cylinder). 

We now look for 0(4) invariant instantons of this theory, i.e. solutions for which 
<p depends only on y. In view of the relations 



tanh y = — cos ( = — yl — r 2 cos r , (25) 

following from fllQ|) , (p3|) , (|i~4"|), and (|15"|), we expect any such solution to have at 
least two turning points with respect to the Euclidean time r, at r = and r = ir. 
Equation of motion for such solutions is obtained from (^) and reads 

- d 2 J + 2 tanh y d y <f> + — -~- f-/i 2 + 3 ) = . (26) 

cosh y K ' 

It is useful to consider also the field <$ = <$/ coshy, for which the equation has no 
"friction" term: 

-<9^ + 0-^^0 + 3 = O. (27) 
cosh y 

For solutions describing tunneling, we require that <fi goes to constant values at y — > 
±oo, and these constant values are different. 

At large /i, eq. ( PB[ ) has a solution for which (f> changes from <p w — to ~ /i in 
a narrow (of order l//x) region near y = 0. This solution is well approximated by the 
kink of double-well quantum mechanics: 

<p w /itanh^/v 7 ^) . (28) 

Because of the finite volume of the de Sitter four-sphere, this solution has a finite 
Euclidean action, of order /z 3 / A. On the other hand, if we formally allow fi 2 to take 
negative values, the usual mechanical analogy, in which is viewed as position of a 
"particle", and j/asa fictitious "time", shows that for /i 2 < —1, eq. (1271) cannot have 
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any bounded solutions except = 0. We conclude that a ^/-dependent solution must 
appear at some critical value /i 2 = // 2 r > — 1. 

This critical (bifurcation) point can be found by the standard argument [EJ, which 
we briefly review here. We assume that at /i 2 just above /i 2 r the new solution is 
confined to a vicinity of = 0. Then, we can linearize eq. (12(1) into 



u 2 

-d 2 + 2tanh?/ <9 V ^-0 = . (29) 

y cosh y 

For /i 2 = 4, this equation has a solution: 

<f) = Atanhy, (30) 



with an arbitrary amplitude A (which should be small to justify using ( P9|) in the 
first place), suggesting that 

/4 = 4- (31) 

A numerical study of the full eq. fl27|) confirms this suggestion: for /i 2 < 4, we have 
found no solutions with a nonzero value of "velocity" d y (j) at y = 0, while for /i 2 > 4 
such solutions exist, with the field approaching (p8|) at large \i. We have also verified 
numerically that the y-dependent solution has a lower action than the Hawking- Moss 
instanton (which in this model is 0o = 0). 

In terms of the static coordinates fll4j) , (|15|), the small-amplitude solution (|30| ) has 
the form 

= -AVl - r 2 cos t , (32) 

corresponding again to small harmonic oscillations near the Hawking-Moss instanton. 
Condition ([H]) is a special case of the condition 



- K"(0o) = 4 , (33) 

obtained in refs. [I], |3| for a wide class of potentials. Here V" is the second derivative 
of the potential, computed on the HM instanton O . The difference with the strongly 
supercooled cases considered in fl], U and in the previous section is that in the present 
case the limit fi ^> 1 leads not to a small, nearly flat-space instanton but to a solution 
like fl28|), for which the "time" duration (in r) is indeed small, of order but the 
spatial extent (in r) is still large, of order of the Hubble radius. 



4 Continuation to real time 

The existence of two turning points, with respect to "time" r, for any 0(4) symmet- 
ric "time" -dependent instantons suggests that these instantons can be analytically 
continued to the Lorentzian signature, and the analytic continuation will describe 
some real-time process. Continuation of time r alone, leaving the radial coordinate 
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intact, is not convenient for studying the real-time evolution, for a reason that will 
be explained shortly. So, we make instead a transformation that acts on both r and 
r, and such that the new time variable is the conventional conformal time rj. (Using 
conformal time for specifying the ^re-tunneling quantum state was advocated in ref. 
00 

The required transformation of coordinates is 

r = -p/ij , (34) 
t = - l -H V 2 -p 2 ). (35) 

where p is a new radial coordinate. The metric takes the form 

ds 2 = i (-drf + dp 2 + p 2 dQ 2 2 ) , (36) 

which is the usual real-time conformal metric of de Sitter space. On the real-time 
segment of the evolution, rj can take any value from — oo to 0. On the Euclidean 
segment, rj is complex (as it should be; note that p is also complex there), but both 
r] and p become real at the turning points r = and r = it. 
At the turning points, eq. ( |3"5] ) gives 

V 2 -P 2 = l, (37) 

which defines a hyperbola in the (rj, p) plane. To each point on this hyperbola, one can 
associate via fl3"l"| ) a unique value of the static coordinate r, satisfying < r < 1, see 
Fig. 2. This means that periodic instantons, such as (0), fl28|), or (|32|) , completely 
determine the field and its derivatives on the hyperbola and thus supply initial data 
for the real-time problem. 

Note that —ir, which would be the time variable if we were to continue r alone, 
goes to infinity at the line p = —rj. It is for this reason that such analytical contin- 
uation is not convenient for studying the real-time evolution. It could still be used, 
though, if we were interested only in the immediate vicinity of the hyperbola fl3"T|). 

It would be very interesting to run the real-time evolution for systems like those 
discussed in the previous sections, with the initial data on the hyperbola determined 
by the turning points that we have found. For our present purposes, however, we 
need only some general features of the evolution, and these can deduced by relatively 
simple arguments. In particular, we will be interested in the general shape of the field 
at a constant rj section, rj = tjq (the vertical line in Fig. 2). 

Indeed, the turning points of (|T8|), fl28|), or (^) do not look anything like droplets 
or bubbles of the new phase. For example, the turning points we have obtained for 
the double-well potential do not interpolate between, say, (f> > at r —> and < 
at r — ► 1. Instead, they all tend to (f> = at r — > 1 (i.e. on the large-p part of the 
hyperbola). We will see, however, that the field acquires the expected form by the 
time it reaches 77 = 770- 
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Figure 2: Hyperbola fl3T|), on which the field is at a turning point, and a section 
of constant conformal time rj. Evolution of the field between these two sections is 
discussed in the text. The line p = —r]r associates a value of the radial coordinate r 
of the static coordinate system to a given point on the hyperbola. 

For definiteness, we take 7] « —1, which is near the top of the hyperbola. Then, 
for p close to zero the field is more or less the same as it was at the turning point, 
because the small p region has been evolving only for a short time. The situation 
is entirely different for large values of p. These regions have evolved for a long time 
before they arrived at rj = r]o, and to find how the field looks there, we have to 
consider not only the turning point itself but also its derivative with respect to rj. 
(This derivative of course does not vanish at the turning point — it is the derivative 
with respect to r that does.) 

For example, the small-amplitude solution ( |32"D in the new coordinates has the 



We need to keep both branches of the square root if want to include both turning 
points, because rf changes its phase by 2n during a half-period of the periodic in- 
stanton. As far as the real-time evolution is concerned, however, we need to consider 
only one of the turning points, and only real rj. Then, it is sufficient to pick only 
the branch corresponding to the turning point in question. In what follows we use 
\fq I = rj < 0. Then, the solution fl38|) in the vicinity of the hyperbola (fffi) is approxi- 
mately <p ~ —A/ 'rj at large \rj\, while its derivative is d v <ft ~ — A. Note that, unlike the 
field, the derivative does not decrease with \rj\ and is of the opposite sign. In other 



form 




(38) 
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words, in the regions with large p the field initially sits close to the top of the barrier 
but is pushed over the barrier, to the "wrong" side. The field in these regions will 
therefore relax to the "wrong" minimum. 

We conclude that the field that was at a turning point on the hyperbola (j37|) 
becomes a droplet of the new phase at r\ = r] , i.e. there is a region with <fi > 
surrounded by <fi < (or vice versa). For the case p ^> H, which is of main interest to 
us in this paper, we expect that the field inside and outside the droplet is close to ±/i. 
The size of the droplet in this case is much larger than the Hubble radius. Indeed, 
the large-/! solution (^) near the part of the hyperbola corresponding to large \rj\ is 
approximately of the form fl38|) but with A replaced by p 2 /^/2. To end up on the 
"wrong" side of the potential (in this case, at <p < 0), the field has to overcome a 
potential barrier of energy density proportional to p 2 a 2 <p 2 ~ p 6 /2i] A , where a is the 
scale factor, while the available kinetic energy is proportional to (9^0) 2 ~ p 4 /2- Thus, 
the field can get to <p < only in regions that started out at \rj\ ~ yfp or larger, i.e. 
the size of the droplet (in p) is of order ^fp/H 3 / 2 . As seen from (p6|) , at i] = — 1 the 
scale factor is equal to 1, so the above estimate should be interpreted as the physical 
size at which the droplet nucleates. 

Since the argument leading to this picture relied on the existence of a nonzero 
derivative with respect to rj of the analytically continued instanton, one cannot help 
wondering if this argument can be adapted to the Hawking-Moss case, when all the 
derivatives are zero. In that case, we do not expect any definite droplet to emerge at 
7] = 7]o, because for p < p CT ~ H (when the HM instanton is relevant) the stochastic 
dynamics ||J] will lead to a random walk of the field at super-Hubble scales. So, 
instead of talking about a specific droplet, we should be talking about the statistics 
of a random field. If we consider initial conditions for which the field is equal to the 
HM instanton on the hyperbola (|37|), but allow for a small nonvanishing "velocity" 
(a derivative with respect to rj), the above argument will go through. Moreover, it 
may be possible to obtain the statistics of the random field at rj = rjo from a suitable 
statistics of the small velocity, much like one obtains a nucleation rate in flat space 
by summing over all different velocities the system may have at the top of the barrier. 
Here, however, we do not pursue this idea any further. 

In the strongly supercooled case of Sect. 2, the droplet will turn into an expanding 
bubble of the new phase.Q We now need to develop the corresponding picture for 
degenerate (or nearly degenerate) potentials. As before, we concentrate on the case 
p 3> H, when transitions between = ±p in a given Hubble volume are due to 
tunneling. (In the case p < p cr ~ H, they can be regarded as overbarrier transitions, 
due to the "Brownian motion" of the field's large-scale component M.) 

In flat spacetime, for a degenerate double- well potential such as (|2ll) , any droplet 
of the phase <t> = p nucleated inside the phase (f> = —p will collapse due to the surface 

* Adding a sufficiently large constant to the potential in (P will allow the field <f> itself to drive 
a (quasi) de Sitter expansion. In that case, the bubble is precisely of the type contemplated in the 
new inflationary scenario jl9|, . 
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tension. We define a nearly degenerate potential by the condition 

R ciit < H- 1 , (39) 

where R CT n is the size of the critical droplet in flat space; if a flat-space droplet is larger 
than the critical size it will expand. In de Sitter space, a droplet with a size much 
smaller than the Hubble radius behaves essentially as its flat-space counterpart. So, in 
either of the above cases, such as a droplet will collapse. We have seen, however, that 
for // ^> H, the size of the droplet is much larger than the Hubble radius. For such 
a droplet, the surface is already causally disconnected from the center. By causality, 
such a droplet cannot collapse. We therefore expect (and will confirm by a numerical 
calculation) that it will develop into an expanding bubble of the new phase, much 
like its counterpart in the supercooled case. (The difference remains that the spatial 
size of the droplet is now strongly dependent on the Hubble radius.) 

The persistence of super-Hubble droplets in the case of a degenerate potential can 
be directly tested by a numerical integration of the real-time equation of motion 

<j> m + - (P y~ p + aV(0) = , (40) 

where a = —l/Hrj is the scale factor, and TC = aH is the comoving Hubble parameter. 
For the model (|2l|) the potential is 

V{<P) = - l -^ + \<P*. (41) 

Fig. 3 presents results of numerical integration of eq. (fl0| ) with the initial condi- 
tions 

<f>(p,V = Vo) = 2exp(-p 2 /0-l, (42) 
<P, v (p,V = Vo) = 0. (43) 

These initial conditions contain a single parameter £, which regulates the spatial 
size of the droplet. Although we have seen that for fi ^> H droplet sizes are quite 
respectably super-Hubble, to prove the point (i.e. the relevance of the Hubble radius, 
as implied by the causality argument), we use for Fig. 3 a mildly super-Hubble value 
I = 2H~ 1 . We see that the comoving size of the configuration rapidly freezes in, 
which means that in the physical space the size grows at the rate of the expansion. 
In contrast, sub-Hubble droplets (not shown in the plot) collapse. 

The role of the de Sitter horizon is emphasized by a comparison of these results 
with the behavior of the same model in a radiation-dominated (RD) universe. In 
the latter case, we expect that even if a droplet starts out as super-Hubble, it will 
eventually become sub-Hubble, as the growing Hubble radius catches up with the 
droplet's size. Once the size of the droplet becomes much smaller than the value of 
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Figure 3: Radial profiles of the field <fi at different moments of conformal time, ob- 
tained by numerical integration of eq. ( 40[) with the potential ([|l|) and initial con- 
ditions (|42|), (f43"D, for p = 100, £ = 2, and r] = —1 (all in units of the Hubble 
rate/radius). Because p is a comoving coordinate, the position of the Hubble radius 
in the plot shifts with time to smaller p, according to ppn = p\r]\- 
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the Hubble radius at that time, the droplet's evolution is essentially the same as in 
flat spacetime. Thus, in an RD universe the droplet will eventually collapse regardless 
of its initial size. These expectations are well born out in numerical integrations of 
eq. © with a = H rj + 1 and H = (77 + l/H )-\ 

The causality argument implies that in a de Sitter spacetime even a droplet of 
false vacuum, in a model with nondegenerate minima, should expand, as long it has 
nucleated with a super-Hubble size. This corresponds to the so-called "jumps up" (in 
potential energy), which are known to occur in de Sitter space both for small fi pl| 
and for /i > fi cr J7J. Our integrations of ([K]), with a nondegenerate V and \i 3> H, 
confirm that the real-time picture implied by causality is correct also in this case, 
although now we find that the expansion of the droplet, at least initially, proceeds at 
a slower rate than the expansion of the universe. 

5 Calculation of decoherence time 

We now have the tools required to understand the influence of a de Sitter horizon 
on quantum coherence, for fields with fi ^> H . If with the potential fl4l|) were 
not a field but a single quantum-mechanical degree of freedom, tunneling would lead 
to quantum oscillations of P(t) = -Pl(^) — Pr{^)i where Pl{r) is the probability to 
observe a negative (positive) 0, i.e. to find the system in the left (right) well. These 
oscillations would be of the form 

P(t) = cos ttt (44) 

with Q proportional to the tunneling amplitude. 

Now that 4> is actually a quantum field in a de Sitter spacetime, the picture is quite 
different. As we have seen, each tunneling event nucleates a super-Hubble droplet of 
the new phase, and the droplet subsequently expands. As a result, the state of the 
field at the center of the droplet changes, say, from <fi « —\h to <ft ~ /■*• Later, another 
tunneling event will occur in the same Hubble volume, with another expanding droplet 
nucleating, and the field changing back to <p ~ — £*■ Although within the given Hubble 
volume the final field becomes indistinguishable from the original, the super-Hubble 
configurations are in most cases quite distinct: if the second droplet has nucleated 
not too soon after the first, there are now two receding bubble walls, while originally 
there were none. If, in such a case, we consider an overlap of the initial quantum 
state with the final one, by taking a product of individual overlaps for all the field 
modes, the rapidly growing difference between the states of super-Hubble modes will 
lead to a rapid decrease in the overlap. One consequence of that will be suppression 
of quantum oscillations of the probability (decoherence). 

To estimate how rapidly the overlap decreases, consider a configuration of two 
bubble walls, with between the walls, and ~ —\l everywhere else. Except for 

the case when the second wall forms very close to the first, we expect the configuration 
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to expand at the rate of the expansion of the universe, so that the comoving volume 
occupied by ~ /x remains to a good accuracy constant, equal to some V. We 
have confirmed this expectation numerically for the special case of two concentric 
bubbles. Let us consider the decrease in the overlap due to the displacement of the 
constant mode of the field from = — /x to = /x in volume V and neglect possible 
displacements of all other modes (radiation). This assumption can only lead to an 
overestimate of the overlap. 

The constant mode has a time-dependent frequency u)q{vi) (with respect to con- 
formal time rj) given by 

u 2 M = 2/xV(77) - 1/v 2 ~ 2/xV(?7) , (45) 

where a = —1/Hrj is the scale factor. Because we consider the case /x ^> H, the 
frequency squared is positive, so there is no amplification of the constant mode, 
and we can simply compute the overlap between the ground states corresponding to 
= ±/x in volume V. To represent the constant mode as a quantized oscillator, we 
switch to the conformal field con f = <pd- The field displacement is now A0 con f = 2/xa. 
With the help of the usual formulas for harmonic oscillators, we find the overlap of 
the two ground states: 

|(0 = /i|0=-/i)| 2 = exp[-Vcu o (A0 conf ) 2 /2] . (46) 

Thus, the overlap decreases to zero as an exponential of a 3 |] 

To obtain the tunneling amplitude, we should explicitly take into account the 
possibility that a bubble can nucleate at different moment of time. Because the 
"static" metric fll6|) is independent of r, there are instanton solutions obtained from 
the previously found instantons by shifting r by a (generally, complex) constant. 
Thus, we consider a shift 

r ^ t - r nucl = r - it nnd , (47) 

with an arbitrary real t nuc i as a parameter. The shifted solution has turning points 
at t = it nuc \ and r = n + it nuc \, so that the turning point hyperbola ( B7p generalizes 
into 

V 2 ~P 2 = exp(-2t nucl ) = • (48) 
The amplitude of tunneling per a small interval of time is now 

QdT nnd = iQdtnacl , (49) 



t This should be contrasted with the behavior of massive gauge theories in fiat spacetime. There, 
a periodic instanton (in theories that have those) describes a transition between two coherent states 
of the field, while the subsequent real-time evolution corresponds to radiation of particles to infinity 
[ p^| . The overlap of the final radiation state with the vacuum, computed for some cases in ref. Q , 
remains finite, although exponentially small, at arbitrarily large times. 
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where Q oc exp(— Se), and Se is the action for half a period of the periodic instanton. 
The quantity 77 nuc i can be interpreted as the moment of nucleation in conformal time. 
Then, then according to (fj8|) , t nuc \ is the same moment in time t = —\n\rj\, which is 
the conventionally defined cosmic time (in units of H^ 1 ). 

Decoherence time can now be calculated by standard techniques. We define Pi,{t) 
and Pr($) as the probabilities to find the field in a chosen Hubble volume to be near 
— /i and [A, respectively. Following ref. || we write P(t) = Pi(t) ~ Pn(t) as 

P(t) = E ("l) n ^ 2n / dt 2n / dt 2n ^ ... dhFfr, t 2 , . . . , t 2n ) . (50) 

n=0 J ° J ° J ° 



Here F(ti,t 2} . . . ,t n ) is the "influence functional" [ Z2| , equal to the overlap between 
two quantum states, obtained from a given initial state (vacuum, in our case) in the 
course of two different bubble nucleation histories. The field in the chosen Hubble 
volume is near ±/x for all intervals between the nucleation times t\,t 2 , . . . , t 2n . How- 
ever, while for t 2i < t < £2*4-1 the field is in the same well for both histories (a diagonal 
state), for t 2 j-i < t <hi the wells are different (an off-diagonal state). 

We have seen that if the system remains in an off-diagonal state long enough 
to allow formation of two well-separated bubbles, corresponding to transitions at 
t = t 2 i_\ and t = t 2 i, the overlap rapidly decreases with t. This means that the main 
contribution to ( ^0] ) comes from pairs of histories for which the transitions at £ 2 i— 1 
and t 2i are "confined" , i.e. occur back to back. For such pairs of histories, we can 
approximate the influence functional as 

n j 

F(h, h, . . . , t 2n ) = n V7<Ht 2j - ta-i) , (51) 
i=i M 

where M~ l is a characteristic timescale, which will be discussed shortly. Eq. ( pT] ) is 
analogous to the "noninteracting blip" approximation of ref. ||. 
Substituting Q5"ip into (0), we obtain 



°° O 2 ™ ft r tr 2n /"*4 / \ 

P(t) = / dt 2n dt 2n „ 2 ... dt 2 = exp(-m 2 /M) . (52) 

r\ IVl JQ JO JO 



00 Q 2n rt r t 2n rU 

E( ' ^ 

n=0 

Eq. (|52]) describes a completely incoherent (i.e. non-oscillatory) relaxation of the 
probability, with decoherence time 

t d = M/n 2 . (53) 

This value sets the timescale on which, due to uncorrelated transitions from <fi w —fi 
to ~ fi and back, the probabilities to observe the system in either of these states 
become roughly equal, Pl ~ Pr « 1/2. 

In the back-to-back events, the bubble walls form so close to each other that they 
annihilate, giving rise to radiation of <p waves. (We have confirmed numerically that 
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this is indeed what happens in the special case of closely spaced concentric bubbles.) 
The overlap of the quantum states, corresponding to a history that contains such an 
event and a history that does not, depends exponentially on the number of radiated 
particles. We therefore expect that the main contribution to P(t) will come from 
those back-to-back events for which there is little or no radiation produced. This 
suggests that a good approximation for the timescale M~ l can be obtained if we 
interpret each Q 2 /M factor in Q52| ) as corresponding to a complete periodic instanton, 
with the field going from one turning point to the other and then back. Then, M can 
be computed from the determinant of small fluctuations near the periodic instanton. 

Although these results were obtained for a strictly degenerate potential, we expect 
them to remain practically unchanged in the nearly degenerate case defined by the 
condition fl39|). Indeed, the main distinguishing feature between these two cases, the 
existence or nonexistence of a large critical droplet in flat space, seems to be totally 
irrelevant for the calculation of decoherence time. 



6 Conclusion 

The existence of turning points for 0(4) symmetric instantons in de Sitter spacetime 
may seem almost obvious (in view of relations (p5])). Nevertheless, it is this property 
that distinguishes them from instantons responsible for the ^-vacuum structure in 
gauge theories and two-dimensional sigma- models. For example, the 0(2) nonlinear 
sigma-model contains a single complex scalar field, which is not affected by conformal 
transformation, and if the metric on the two-sphere is chosen in the form (|2]) (where 
now fi, v — 1, 2) the single- and multi-instanton solutions have exactly the same 
form as in the flat Euclidean space. For the latter case, they are given in ref. [[23 



In particular, the one-instanton solution centered at the origin can be written as 
w\ = bz, where z = x 1 + ix 2 , and b is a complex parameter (its modulus determines 
the instanton size). Writing z = e l<?i tan|, where 9 and <p are t ne usual spherical 
angles, we can define the "static" coordinates, analogous to (|I~4"D , ([15]): 

X = cos 6 , (54) 
r = <f>, (55) 

with X running from —1 to 1, and r from to 2n. In these coordinates, the one- 
instanton solution takes the form 

/ 1 Y"\ 1/2 

Wl = be iT . (56) 



1 + X. 

The real and imaginary parts of W\ never vanish simultaneously; hence, the solution 



has no turning points. Similar considerations apply to the BPST instanton f24 
four- dimensional massless gauge theories. 



in 
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Classically, the instanton action, in either of these cases, is exactly the same as 
in the flat space. We expect that the inclusion of quantum corrections will result in 
a tunneling amplitude that is somewhat reduced compared to the flat space. This 
is because, as in the model of Sect. 2, the maximal instanton size is limited by the 
Hubble radius. Thus, for example, for the BPST instanton, we expect the amplitude 
to be suppressed by exp[— 871 2 /g 2 (H)}, similarly to the suppression of instantons at 
finite temperature p5[] . Nevertheless, the resulting ^-dependence of the partition sum 
is nonzero and represents an effect of quantum coherence. In contrast, instantons 
with turning points (periodic instantons) describe, as we have seen above, completely 
incoherent tunneling. 

The existence of turning points suggests that the instanton solutions found in 
Sects. 2 and 3 can be continued to the Lorentzian signature, and the turning points 
will provide initial conditions for the real-time evolution. To explore this real-time 
evolution we have made a simultaneous transformation, eqs. Q3"l|), (|33|), of the tem- 
poral and radial coordinates, such that the new coordinates are the conventional 
conformal coordinates of de Sitter spacetime. In these coordinates, the initial condi- 
tions are imposed on the hyperbola in the time-radius plane, see eq. (|37|) and Fig. 2, 
and our main results were obtained by studying this initial value problem. 

These results are as follows, (i) We have found that the field that starts out as 
a turning point on the hyperbola (|37|) develops, by the time it reaches a section of 
constant conformal time 1] = r] , into a super- Hubble droplet of the new phase, (ii) 
We have argued from causality and confirmed by numerical integration that a super- 
Hubble droplet will become an expanding bubble even in cases (such as degenerate or 
nearly degenerate potentials) when there are no expanding bubbles in the flat space 
or such bubbles must have very large sizes, (iii) We have shown that super-Hubble 
effects caused by the expansion of de Sitter bubbles completely wash out any coherent 
quantum oscillations that one might associate with a double-well potential, and we 
have calculated the corresponding decoherence time. 

As we have mentioned in the introduction, these results find an immediate ap- 
plication in cosmological scenarios in which there are many "worlds" that might be 
connected by quantum coherences. What we show, in a nutshell, is that a de Sitter 
(or quasi-de Sitter: for example, inflationary) stage, while enhancing tunneling in 
many cases, at the same time completely destroys any traces of such coherences. In 
this way, the assumption that a de Sitter stage existed in the distant past removes the 
problems of interpretation presented by a coherent superposition of states of the en- 
tire universe. It also alleviates possible worries about catastrophic events that might 
be thought to result if the universe today somehow picked up the "wrong" component 
of a superposition. 

The author thanks V. Rubakov for comments on the manuscript and M. Shaposh- 
nikov for discussion of the results. He also thanks the Institute of Theoretical Physics, 
University of Lausanne, where this work was completed, for hospitality. This work 
was supported in part by the U.S. Department of Energy through Grant DE-FG02- 



19 



91ER40681 (Task B). 



References 

[1] S. W. Hawking and I. G. Moss, Phys. Lett. HOB, 35 (1982). 

[2] A. H. Guth and E. J. Weinberg, Nucl. Phys. B212, 321 (1983). 

[3] L. G. Jensen and P. J. Steinhardt, Nucl. Phys. B237, 176 (1984). 

[4] A. A. Starobinsky, in: Field Theory, Quantum Gravity, and Strings, eds. H. J. 
de Vega and N. Sanchez (Springer, 1986). 

[5] A. S. Goncharov and A. D. Linde, Sov. J. Part. Nucl. 17, 369 (1986). 

[6] N. Deruelle, Mod. Phys. Lett. A 4, 1297 (1989). 

[7] V. A. Rubakov and S. M. Sibiryakov, Theor. Math. Phys. 120, 1194 (1999) fgr] 

| qc/990509311 . 

[8] S. Gratton and N. Turok, Phys. Rev. D 63, 123514 (2001). 

[9] S. Chakravarty and A. J. Leggett, Phys. Rev. Lett. 52, 5 (1984); A. J. Leggett 
et al, Rev. Mod. Phys. 59, 1 (1987). 

[10] S. Habib, Phys. Rev. D 46, 2408 (1992). 

[11] C. G. Callan, R. F. Dashen, and D. J. Gross, Phys. Lett. 63B, 334 (1976); R. 
Jackiw and C. Rebbi, Phys. Rev. Lett. 37, 172 (1976). 

[12] A. D. Linde, Phys. Lett. B 175, 395 (1986). 

[13] G. W. Gibbons and S. W. Hawking, Phys. Rev. D 15, 2752 (1977). 

[14] S. Coleman and F. De Luccia, Phys. Rev. D 21, 3305 (1980). 

[15] S. Fubini, Nuovo Cim. 34A, 521 (1976); L. N. Lipatov, Sov. Phys. JETP 45, 216 
(1977). 

[16] S. Khlebnikov, V. Rubakov, and P. Tinyakov, Nucl. Phys. B367, 334 (1991). 
[17] F. R. Klinkhamer and N. S. Manton, Phys. Rev. D 30, 2212 (1984). 
[18] G. 't Hooft, Phys. Rev. D 14, 3432 (1976). 
[19] A. D. Linde, Phys. Lett. 108B, 389 (1982). 

[20] A. Albrecht and P. J. Steinhardt, Phys. Rev. Lett. 48, 1220 (1982). 

20 



[21] A. D. Linde, Phys. Lett. B 131, 330 (1983). 

[22] R. P. Feynman and F. L. Vernon, Ann. Phys. (N.Y.) 24, 118 (1963). 

[23] A. A. Belavin and A. M. Polyakov, JETP Lett. 22, 245 (1975). 

[24] A. A. Belavin et al, Phys. Lett. 59B, 85 (1975). 

[25] D. J. Gross, R. D. Pisarsky, and L. G. Yaffe, Rev. Mod. Phys. 53, 43 



21 



